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We point out that naive pentaquark correlations function include two-hadron-reducible contribu- 
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tions, which are given by convolution of baryon and meson correlation functions and have nothing 
to do with pentaquark. We show that the two-hadron-reducible contributions are large in the 
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operator product expansion of the correlation functions of three existing works on the pentaquark. 

o ' 

■«Hj- , Therefore, it is dangerous to draw a conclusion from the sum rules using naive pentaquark correla- 



tion functions with naive ansatz for the spectral function under the dispersion integral. Instead, we 
propose to use the two-hadron-irreducible correlation function, which is obtained by subtracting 
the two-hadron-reducible contribution from the naive correlation function. Taking one of the works 
as an example we demonstrate how drastically the results can change if we remove the two-hadron- 

•rH , 

reducible part from the naive correlation function. We obtain the result opposite to the original 
work for the parity of the pentaquark. 
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Possible existence of an S = +1 exotic baryon state has recently been reported by LEPS 
collaboration in Spring-8 In this experiment, a peak with the mass ~ 1540 MeV and 
the width bounded by 25 MeV was observed in K + n channel from 'jn — > K + K~n. It was 
confirmed by subsequent experiments □□□□□□□□□ . In all cases the mass is near 
1540 MeV and the width is small. 

This state cannot be a three-quark state since it has S = +1, and the minimal quark 
content is (uudds). It has come to be called "pentaquark + ". Other quantum numbers, 
spin, isospin and parity, have not yet been experimentally determined. Besides the calcu- 
lation in the chiral quark soliton model ^l], which motivated the pentaquark search by 
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15. 16. 17, 181 have been 



LEPS collaboration, various theoretical approaches 
performed, but no consensus has yet been seen on these quantum numbers. 

Clearly, the discoveries of exotic hadrons have opened a new field of strong interaction 
physics and will lead us to a deeper understanding of QCD. In particular, what role con- 
finement plays in exotic hadrons is an extremely interesting and important issue. 

In this paper, we focus on the application of the QCD sum rule approach to the 
pentaquark and discuss an issue which is characteristic for exotic hadrons. In QCD sum 
rule approach, the correlation function of interpolating fields is calculated by the use of 
the operator product expansion (OPE), and is compared with the spectral representation 
via dispersion relation. The sum rules relate hadron properties to the vacuum expectation 
values of QCD operators, such as (0|gg|0) and (0|(o¡ s /7r)(j 2 |0). 

Uü to now, three groups have reported results of QCD sum rules for the pentaquark 
0, 0, Q. zhu Q estimated the mass of the pentaquark states with QCD sum rules 
and found that pentaquark states with isospin I = 0,1,2 lie close to each other around 
1.55 ± 0.15GeV. He used the interpolating fields for the I = 0, 1, 2 pentaquark states where 
Lree quarks and the remaining sq are both in a color adjoint representation. Matheus et al. 
"y\ used a linear combination of two independent interpolating fields: one is made of two 
scalar ud diquarks and the other of two pseudo-scalar ud diquarks, together with s. Their 
result for the pentaquark mass is me+ = 1.55 ±0.10 GeV. It should be noted, however, that 
their interpolating field has isospin 1 = 1 although it was stated to have isospin / = in 
Ref. 17 1 . Sugiyama, Doi and Oka Jj| employed an interpolating field with J = 1/2, 1 = 
and S = ±1 constructed from color anti-triplet scalar and pseudoscalar ud diquarks and an 
s quark. They derived sum rules for the positive and negative parity states. They found 
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that the pole residue obtained from the sum rule is positive for the negative parity state but 
negative for the positive parity state, which implies that the obtained negative parity state 
is a real one but the pole in the positive parity spectral function is spurious. Thus, they 
concluded that the parity of the pentaquark is negative. The mass me+ was estimated to 
be near m e + ~ 1.5 GeV. 

Although these three groups obtained the pentaquark mass close to the observed value, 
we point out that their naive pentaquark correlation functions include two-hadron-reducible 
contributions, which are due to noninteracting propagation of the three-quark baryon and 
the meson and therefore have nothing to do with the pentaquark. These contributions exist 
in the correlation function only for exotic hadrons and are potentially large. Therefore, it 
is dangerous to draw a conclusion from the sum rules using naive pentaquark correlation 
functions with naive ansatz for the spectral function under the dispersion integral. Instead, 
we propose to use the two-hadron-irreducible correlation function, which is obtained by 
subtracting the two-hadron-reducible contribution from the naive correlation function. We 
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show how large the two-hadron-reducible contributions are in t 
functions considered in the above three works. Taking Ref . jl8 1 
demonstrate how drastically the results change if we remove the two-hadron-reducible part 
from the naive correlation function. 

In order to derive QCD sum rules for the pentaquark, one considers the correlation 
function 

Upip) = -if d 4 xe^(T(77p(x)77p(0))), (1) 

where r]p(x) is the interpolating field for the pentaquark, a composite operator made of five 
quark fields. The spectral function is given by the imaginary part of the correlation function 
and represents the physical spectrum generated by the interpolating field as 

Pp(p) = ^5( P - Pn )(0\r)p(0)\n)(n\f) P (0)\0) (p° > 0). (2) 

n 

In the QCD sum rule approach the spectral function is usually parametrized by a pole plus 
continuum contribution, which is adopted also for the pentaquark in Refs. 
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p P (p) = Z5(p 2 - M 2 P ) + 6(p 2 - s th )p° P ^(p 2 ), (3) 

where Mp is the mass of the pentaquark, s t h the effective continuum threshold energy 
squared and Pp PE (p 2 ) the imaginary part of the correlation function in the OPE. For or- 



dinary hadrons this procedure is probably all right but for exotic hadrons it might be too 
naive as we discuss. 

A remarkable feature of the pentaquark is that it can be decomposed into a color-singlet 
three-quark state, baryon, and a color-singlet quark-antiquark state, meson. (Hereafter, we 
use the term baryon when its minimal quark-content is qqq.) Therefore, the interpolating 
field for the pentaquark can be expressed as a sum of the product of baryon and meson 
interpolating fields: 

Vp(x) = J2vU x )Vm( x ) (4) 

i 

where rf B (x) and w M (x) are color-singlet baryon and meson interpolating fields, respectively: 

Vb ~ £ab C q a q b q c , vm ~ q e q e - 

Due to this separability, the pentaquark correlation function, Eq.ífTj). has a part in which 
the baryon and the meson propagate independently without interacting each other. We 
define this part as the two-hadron-reducible (2HR) part and the rest of the correlation 
function as the two-hadron-irreducible (2HI) part. They are respectively given by 

(T( VP (x)fjp(ow HR = E(nvB^)ém}(nvM^)vM(o))), 

(T( VP (x)fj P (0))) 2HI = ^(^(^^(O)))-^^^^)^^)))^^^)^^))) (5) 

ij 

and 

n 2 P HR (p) = ¿E/ ^m(q)U%(p~q), 

U 2 p HI (p) = U P (p)-ij: J ^n%(q)U%(p-q), (6) 

where U l ¿(p) and U M (p) are baryon and meson correlation functions: 

U%(p) = -ij éxév%T(rf B (x)ff B (0))), (7) 
lV M (p) = -iJ d^(T(rf M (x)rf M (0))). (8) 

Diagrammatically, the 2HR and 2HI parts of the correlation function are represented as 
Figures ^ and El respectively. Clearly, the 2HR part of the pentaquark correlation function 
is completely determined by the baryon and meson correlation functions and has nothing to 
with the pentaquark. 
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FIG. 1: Two-Hadron-Reducible (2HR) diagram. 




Some comments are in order here. The 2HR and 2HI parts have the same analytic 
property as the total correlation function. The concept of the two-hadron-reducible (2HR) 
and two-hadron-irreducible (2HI) contributions in the correlation function of the composite 
particle is a generalization of the two-particle-reducible (2PR) and two-particle-irreducible 
(2PI) contributions in the correlation function of the elementary particle. 

Let us next look at the separation of the 2HR and 2HI parts in the spectral function. We 
suppose that the lowest states generated by r¡ B and r¡ M are spin-1/2 baryon B and spin-0 
meson M, respectively. Consider only the contribution of the BM scattering states in the 
spectral function just for simplicity, which is given by 

p B P M {p) = E / d 3 kd 3 q{2n) 3 5\p -k- q)(0\r] P (Q)\kqsout)(kqsout\f¡ P (0)\0) (p° > 0),(9) 

s J 

where q and k are the momentum of the baryon and the meson, respectively, and s the 
spin projection of the baryon. \kqsout) denote the baryon-meson scattering states with out 
boundary condition. One can relate the matrix elements in Eq.(|ü¡) with the T-matrix for 
the BM scattering by means of the reduction formula as 

(kqsout\fj{0)\0) = X b Xm 1 , 

í f d^q' 1 1 1 

xü{qs) ll — i —^T(k, q,k + q-q', q')j7- ^ j-r -7. „ . \ , (10) 

[ J (27r) 4 [k + q — q'Y — m ¿ + ir¡ § — M + ir¡ J 
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(0\r¡(0)\kqs(mt) = \* B \* M ¡ 



x { 1 + i I -FTÜJn ^ 2 ■ 7, It ■ Tt ( fc ' + 9 - V', </)} u(qs)(ll) 

(27TJ 4 (k + q — q') z — m ¿ — ir¡ $ — M — ir¡ I 



where Uk = -y/k 2 + m 2 , E q = ^/q 2 + M 2 , \b(m) is the coupling strength of the interpolating 
field to the baryon (meson) state and M(m) is the baryon (meson) mass. T(fc, g, fc', q') is 
the T-matrix for the process kq — >■ fc'g'. By substituting Eqs. ífTÜj) and (fTTj) into Eq.Q, we 
obtain 



(i 4 g 1 1 



(27r) 4 (p — q) 2 — m 2 + ir¡ $ — M + ir¡ 

d^g <¿y 1 1 

(2vr) 4 (2vr) 4 (p - g) 2 - m 2 + ir¡ $ - M + ir) 

xT(p-q,q,p- q',q')j — — — — — \ . (12) 

f — M + irj (p — q J 2 — m ¿ + ir¡ ) 

The first term of Eq. (jl2|) is the 2HR contribution due to the trivial noninteracting contri- 
bution of the BM intermediate states and the second term is the 2HI contribution. If the 
pentaquark is a bound state in the BM channel the spectral function has an additional con- 
tribution from the bound state, while if the pentaquark is a resonance in the BM channel 
the effect of the pentaquark lies in the BM T-matrix as a pole at a complex pentaquark 
energy. In any case, the 2HR contribution is not related to the pentaquark. 

Some further comments are in order here. 2HR contributions discussed here exist com- 
monly in the correlation functions for exotic hadrons but not for ordinary hadrons. Crucial 
assumption here is confinement. Namely, we assume that only color-singlet states contribute 
to the spectral function. Therefore, the separability of the pentaquark into color-singlet 
baryon and meson is the origin of the existence of the 2HR contribution. Another point to 
be mentioned is that by removing the 2HR part we do not take out all the baryon-meson 
continuum contributions. We just take out noninteracting contributions. The continuum 
contributions included in the 2HI part cannot be separated from the pentaquark contribu- 
tion. The 2HI part of the spectral function does not have positivity because of subtraction, 
but the positivity is not necessary to construct sum rules. 

Let us turn to the separation of the 2HR and 2HI parts in the OPE. There are various 
ways of constructing the interpolating field for the pentaquark. So far, three ways of the 
construction have been proposed. Zhu used the following interpolating field for the 



6 



/ = pentaquark: 

1 



%hu = ^e abc {u Ta C l5 d b )[u%sH l5 d c ) -{u^ d)\, (13) 



where a, b, c, . . . are color indices and C = «7270- Mathues et al. |Tj| and Sugiyama et al. 
considered the interpolating fields based on the conjecture by Jaffe and Wilczek 

^Matheus eí al. = + ??2, 

771 = -^e afec (u Ta C , 7 5^)(n Tc C'75á e )C'(s e ) T -(n^d), 

r]2 = -^e abc {u Ta Cd b ){u Tc C l5 d c )C{s a f -{u++ d), (14) 

Í/Sugiyama eí al = € ahc € de f€ cfg {u Ta Cd!') {u Td C l5 d C )C {s 9 ) T . (15) 

We first note that after Fierz transformation the interpolating field can be expressed as a 
sum of the product of baryon and meson interpolating fields. The interpolating field used 
by Zhu can be written as 

Vzhu = -^e abc {u Ta C l5 d b ) {-d c {si l5 u) - ^d^si^^u) 

-\a^d c {sia^ l5 u) - Y l5 d c {si lfl u) - l5 d c {siu) - {u <-> d)\ . (16) 

The interpolating fields considered by Matheus et al. and Sugiyama et al. are also written 
as follows, respectively: 



í?Matheus eí al. — tf]i + T] 2 

r¡i = -^e abc {u Ta C l5 d b ){u Tc C l5 d c )C{s c ) T -{u^d) 

= ~ {-[e abc {u Ta C l5 d b )u c }{s l5 d) - [€ abc {u Ta C l5 d b )^u c ]{s l5llí d) 

+ \[e abc {u Ta C l5 d b )a^u c }{sa flvl5 d) + [e^u^C^^u^s^d) 
-[e abc {u Ta C l5 d h ) l5 u c ]{sd)} -{u^d) 
V2 = -^e abc {u Ta Cd b ){u Tc Cd c )C{s 9 ) T - {u ^ d) 

= ~ {-[e abc {u Ta Cd b )u c }{sd) + [e abc {u Ta Cd b )^u c }{s lfi d) 

+ \[e abc {u Ta Cd b )a^u c }{sa fllJ d) + [e abc {u Ta Cd b )Y l5 u c }{s l5lfl d) 
-[e abc {u Ta Cd b ) l5 u c }{s l5 d)} -{u^ d), (17) 

^Sugiyama et al. = ^€ abc {u Ta C d°) {~d C {s l5 u) - ^ d° {s l5l ^u) 

+ \a^d c {sa fiul5 u) + 7^75^(57^) - l5 d c {su) - {u c¿)} . (18) 



The 2HR parts in the correlation functions are obtained by contracting a quark field in 
rf B (x) with that in ff B (0) and contracting a quark field in rf M (x) with that in rf M (0), while 
the terms including contraction of a quark field in the baryon with that in the meson belong 
to the 2HI part. We show the total and the 2HI part of the spectral function in the OPE. 
For the correlation functions considered by Zhu, the results are 

í (V 2 f 11 (P 2 ) 3 , v (V 2 f , v 

Pzhu(p) = { 2 2o 3 , 5 , 7?r s ~ WZW ms{qq) + W ms{SS) 

2 19 3 2 57r 6 ^7r 1 2 n 3 2 7r 4Wyyy/ 2 12 3% 4 wy 7 j F 

+ 3203252-8 m * 2 16 3 2 5vr 6 m) 2 17 3 2 5vr 6 ^ ; 
25(p 2 ) 3 7(p 2 ) 3 

1 (/)5 ¡™. + J^<*_>-J^<*> 



2 20 3 2 5 2 tt 8 4 2 17 3 2 5tt 6 2 17 3 2 5vr 6 
5Qo 2 ) 3 7(p 2 ) 3 

~~~~~~ {q9sCr ' Gq) ~ ~~~~~ {s9s(J ' Gs) ' (19) 

where we have included the condensate (^G 2 ) , which is missing in Ref . }l6_ . For the 
correlation functions considered by Matheus et al. and Sugiyama et al. the results are 

n _í (p 2 ) 5 , (P 2 ) 3 (P 2 f 

PMatheus eí ol IPJ - | C l 21832527-8 + 2 14 3 2 5vT 6 ^ 2 17 3 2 5tT 6 ^ V ' 

(-2)2 _ ^2)2 _ 1 

~ Cl 2~~~~ ms{s9(r ' Gs) + C3 2^^ {qqqq) j^ 

+Cl ~~~W ms ~ C1 2W (SS> + C1 2W (SF ' GS) ' 
2Hi / x í ÍP 2 ) 5 (p 2 ) 3 , . (p 2 f ,a a ^ 

P M athcu S et ol (P) ~ IC!^^^^^ + C l 2 l 63 2 57r 6 m ^ SS ) ^ffl^V^ } 

(P 2 ) 2 - (P 2 ) 2 - 1 

" Cl 2W ms(sF ' Gs) + c *2~~j~~ {qqqq) ) lt> 

(p 2 f (p 2 f _ (p 2 f 

+Cl 2W ms " C1 2W (SS> + C1 2W (SF ' GS) ' (20) 

where c_ = 5t 2 + 2í + 5, c 2 = (í - l) 2 , c 3 = 7í 2 - 2t - 5 and c 4 = llt 2 + 14í + 11, 

f (p 2 ) 5 (p 2 ) 3 Mj _W__A r , 

PSugiyama ei ol IPJ j 2 16 3 2 5 2 7tT 8 2 12 3 2 5tT 6 ' 2 14 3 2 5tT 6 ^ 7T ' 

¡m s (sg s a ■ Gs) \$ 



2 13 3% 6 

/ p 2)5 / p 2N4 _ _p2)3 

+ 2^3W ms ~ 2W (SS> + 2W ' GS) ' 



(V) = f {p2)5 - 7{p2)3 m (ss) + {p2)3 < a * 

3 ugiyamaeUÍ W 2 19 3 2 5% 8 2 15 3 2 5vr 6 " S 7 2 18 3 2 5tT 6 V 7T 



7p 4 _ } 

7(p 2 ) 5 7(p 2 ) 4 7(p 2 ) 3 

l™* + ^^{ss)-^-Asg s a-Gs). (21) 



2 19 3 2 5 2 tt 8 2 16 3 2 5tt 6N ' 2 16 3 2 tt 6 

As can be seen, the 2HR part of the correlation function is large in general at least of the 
same order as the 2HI part. In particular, for the interpolating field, Eq. (jl5|) . the Wilson 
coeffcients of the operators in the 2HR part are —15/7 of the 2HI part up to dimension 6 
except for the operator, ^G 2 . 

Now, taking the interpolating field of Sugiyama et al. ¡181] . Eq.(|lñ¡). we will demonstrate 
how drastically the results of the sum rule can change if we remove the 2HR part from 
the total correlation function. First, we explain the results of Ref . [ísl] . They defined the 
positive-parity and negative-parity spectral functions by 

P±(Po) = -Tr[{ lü + l)p( Po ,p = 0)]. (22) 

By parametrizing the spectral function by a pole plus continuum contribution as 

P±(po) = \X±\ 2 6(p - m±) + 6(p - ^)p°p PE (Po), (23) 



they derived sum rules, 

iA±|2exp (-¿)=r^ FEwexp (-±)- <24) 

|A±|2m2±6XP (-it) = jT^'^exp , (25) 

where ^/s~th is the effective continuum threshold energy and the parameter of the weight 
function. Meorei is called the Borel mass. FigOl shows the right-hand side of Eq. tfl^ as 
a function of MBorei with the standard values of the quark mass and QCD condensates, 
m s = 0.12 GeV, (ss) = -0.8 x (0.23 GeV) 3 , (fG 2 ) = (0.33 GeV) 4 , m 2 = 0.8 GeV 2 , where 



m 2 , = (qg s (T ■ Gq)/(qq) and y/s~¡j¡ is taken to be 1.8 GeV. The solid and dashed lines are for 
the positive and negative-parity states, respectively. As can be seen, the right-hand side of 
Eq.(|21J) is positive for the negative-parity state but negative for the positive-parity state. 
Since the left-hand side of Eq.(¡24p must be positive, it was concluded that the obtained 
negative parity state is a real one but the pole in the positive parity spectral function is 
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spurious. The mass of the B + was also predicted to be ~ 1.5 GeV from the ratio of Eq. (j25|) 
and Eq.(El. 

n 

Though we have not been fully convinced if the procedure of Ref.[18| is legitimate, we 
temporarily accept it in order to focus on the effect of removing the 2HR contributions 
on the results of sum rules. We replace the total spectral function by the 2HI part but 
assume the same form, Eq.(|23|). for the spectral function under the integral, where pp PE 
in the continuum contribution is also replaced by the 2HI part. FigE] is the same figure 
as Fig|31 after replacing the total spectral function by the 2HI part. Now, the right-hand 
side of Eq.(J2I|) is positive for the positive-parity state but negative for the negative-parity 
state. This result was expected because the Wilson coeffcients of the operators in the 2HR 
part are —15/7 of those in the 2HI part up to dimension 6 except for the operator, ^G 2 . 
Therefore, the sum rule for the 2HI part of the spectral function leads us to the opposite 
conclusion that the obtained positive parity state is a real one but the pole in the negative 
parity spectral function is spurious. The mass of the B + is obtained to be ~ 1.4 GeV. This 
result shows that it is crucial to remove the 2HR part from the total correlation function 
in order to correctly extract information on the pentaquark. We, however, would like to 
suspend the conclusion that the parity of the pentaquark is positive until we make sure that 
all the procedures of deriving sum rules are reasonable. We will report the results of our 
reinvestigation in near future ¡2c| . 

Some final comments are in order here. Logically, there is nothing wrong to use the total 
correlation function. In order to derive sum rules, however, if the correlation function has 
a large background contribution which has nothing to do with what one wants to extract 
information about, pentaquark in the present case, it would be extremely difficult unless 
one knows very precise information about the background. It is much better if the back- 
ground can be exactly separated, which is what we proposed in the present paper. Also, 
we would like to make a brief comment on the Lattice study of the pentaquark. In the 
lattice calculation, the mass of the pentaquark is extracted from the long-time behavior of 
the total correlation function, which contains large two-hadron-reducible contribution. It is 
necessary to remove this two-hadron-reducible contribution in order to obtain information 
on the pentaquark baryon[15]. It would be interesting to see the relevance of the present 
work in the context of the Lattice study of the pentaquark. 
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FIG. 3: |A| 2 exp(-Mp/M]3 orcl ) vs. M Bore i for the total spectral function with ^/s¡¿ = 1.8 GeV 
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